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ABSTRACT 

In Neurospora crassa, the methionine synthase gene 
met-8 plays a key role in methionine synthesis. In 
this study, we found that MET-8 protein levels were 
compromised in several mutants defective in proper 
heterochromatin formation. Bioinformatics analysis 
revealed a 50-kb AT-rich region adjacent to the met-8 
promoter. ChIP assays confirmed that trimethylated 
H3K9 was enriched in this region, indicating that 
heterochromatin may form upstream of met-8- In an 
H3/Cd/7 mutant strain, the output of mef-Swas dramat- 
ically reduced, similar to what we observed in mu- 
tant strains that had defective heterochromatin for- 
mation. Furthermore, the production of ectopically 
expressed met-8 at the his-3 locus in the absence 
of a normal heterochromatin environment was inef- 
ficient, whereas ectopic expression of mef-S down- 
stream of two other heterochromatin domains was ef- 
ficient. In addition, our data show that the expression 
of mig-6 was also controlled by an upstream 4.2-kb 
AT-rich region similar to that of the met-8 gene, and 
we demonstrate that the AT-rich regions adjacent to 
met-8 or mig-6 are required for their peak expression. 
Our study indicates that mef-Sand mig-6 may repre- 
sent a novel type of gene, whose expression relies 
on the proper formation of a nearby heterochromatin 
region. 

INTRODUCTION 

Large regions in higher eukaryotic genomes are composed 
of repetitive deoxyribonucleic acid (DNA) elements, which 
are preferential targets for heterochromatin assembly. In 



multicellular organisms, constitutive heterochromatin is lo- 
cated at pericentromeric and telomeric regions. An addi- 
tional form of facultative heterochromatin, found in the in- 
active mammalian X chromosome, is characterized by the 
presence of trimethylated H3K27 (1,2). Heterochromatic 
regions are highly condensed in structure and are charac- 
terized by their transcriptionally repressed state, as exem- 
plified by the phenomenon of position-effect variegation in 
Drosophila melanogaster, in which the expression of a eu- 
chromatic gene is compromised when it is relocated near 
or within a block of heterochromatin (3). In animals and 
fungi, histone H3 lysine 9 trimethylation (H3K9me3) is 
a distinct feature of heterochromatic loci, whereas in the 
plant Arabidopsis thaliana, histone H3 monomethylation 
and dimethylation at lysine 9 are the predominant mark- 
ers for silenced heterochromatin. Heterochromatin protein 
1 (HPl) recognizes and binds to H3K9 methylation through 
its chromodomain for heterochromatin assembly in most 
well-studied model systems (4). 

In fission yeast, the protein Rikl interacts with Cul4 
to target the histone methyltransferase Clr4, forming the 
ClrC complex for H3K9 methylation (5-7). In plants and 
metazoans, DDBl is the only Rikl-like protein (8,9) and 
is a component of Cul4-based E3 ubiquitin ligases in eu- 
karyotes (10-12). Knockdown of the ddbl or cul4 genes in 
mammalian cells decreases trimethylation of H3K4, H3K9 
and H3K27 (13). Moreover, DDBl and Cul4 proteins ap- 
pear to associate with trimethylated H3K9 and H3K27 
(13), the principal markers of heterochromatic domains 
and transcriptionally silenced loci. These findings suggest 
that DDBl and Cul4 may be important for gene silenc- 
ing and heterochromatin formation in mammals, in ad- 
dition to their role in the DNA damage repair pathway. 
In Neurospora crassa, the histone methyltransferase DIM- 
5 is required for H3K9 methylation and it functions with 
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HPl in DNA methylation and heterochromatin forma- 
tion at most regions (14-16). Recent studies demonstrated 
that in N crassa, DDBl and Cul4 interact with DIM-5 
via a bridge comprising DCAF26 (DDB1-Cul4- Associated 
Factor 26) and DIM-7 (Defective In Methylation 7). These 
proteins form a complex that is essential for histone H3K9 
trimethylation and DNA methylation at AT-rich DNA re- 
gions, where heterochromatin may potentially form (17- 
20). These results suggest that the Cul4-based histone 
methyltransferase complex may be involved in the regula- 
tion of heterochromatin formation and gene expression in 
N crassa and other higher eukaryotes. 

Many genes from diverse organisms are known to be si- 
lenced by the nearby heterochromatin (3), but the expres- 
sion of certain genes is known to be dependent on hete- 
rochromatin (21-24). In D. melanogaster, these heterochro- 
matic genes are immune to the silencing effects of hete- 
rochromatin and are known to depend on heterochromatic 
locations and SU(VAR) proteins for normal expression (25- 
27). The notion that repetitive elements can create hete- 
rochromatic repressive environments by recruiting silenc- 
ing components is well accepted, but it is unknown how 
this environment might facilitate active expression of het- 
erochromatic genes. The expression of some essential hete- 
rochromatic genes within these repressive regions presents 
a paradox, and the understanding of which could provide 
key insights into the effects of heterochromatin structure on 
gene expression (24). 

To further understand the mechanism of heterochro- 
matin for gene regulation in eukaryotes, we sought to 
identify pathways that selectively regulate the expression 
of heterochromatic genes. Here we report that normal, 
high-efficiency expression of the methionine synthase gene 
(met-8) requires the presence of a proximal 50-kb AT-rich 
region of heterochromatin in N crassa. met-8 expression 
appears to be compromised when proper heterochromatin 
formation was bothered, which is in contrast with the 
well-accepted notion that repetitive elements create hete- 
rochromatic repressive environments by recruiting the si- 
lencing components for gene silencing (28,29). In addition, 
we found that the expression of menadione-induced gene-6 
(mig-6) was positively controlled by its upstream 4.2-kb AT- 
rich region, in a similar way to the met-8 gene. Our findings 
not only strongly implicate Cul4-based complex activity in 
gene expression via its effects on the proper formation of 
heterochromatin but also provide evidence of a clear func- 
tion for heterochromatin in maintaining expression of the 
met-8 and mig-6 genes. 



MATERIALS AND METHODS 

Strains and culture conditions 

In this study, N. crassa 87-3 (bd, a) (30) was used as the wild- 
type strain. The H3K9R knock-in mutant was newly gen- 
erated from the ku70^^ {bd, a) genetic background strain. 
cul4^^, ddbl^^, dcaf26^^, dim-7^^, dim-5^^, hpo^^ and 
dim-2^^ strains, generated previously (17,18), were also in- 
cluded in this study. The 301-6 (bd, his-3, A) and dim-7^^ 
(bd, his- 3) strains were used as the host strains for his- 3 tar- 
geting constructs. 



Liquid cultures were grown at 25°C with shaking in min- 
imal medium (1 x Vogel's and 2% glucose) for ~22 h in con- 
stant hght (LL). 

Protein extraction for 2D gel 

Tissues were harvested by filtration and ground in liquid 
nitrogen. The fine powder was precipitated overnight at 
— 20°C with 10% (w/v) trichloroacetic acid in cold ace- 
tone containing 0.07% (v/v) 2-mercaptoethanol. The mix- 
ture was centrifuged at 40 000 g at 4°C for 1 h, and the pre- 
cipitates were washed with cold acetone containing 0.07% 
(v/v) 2-mercaptoethanol, 1-mM phenylmethylsulfonyl fluo- 
ride (PMSF), l-(jLg/ml pepstatin A and l-jxg/ml leupeptin. 
Dried pellets were then dissolved in a 7-M urea solution 
with 2-M thiourea, 20-mM dithiothreitol, 1-mM PMSF, 1- 
|jLg/ml pepstatin A and l-fxg/ml leupeptin before centrifu- 
gation at 100 000 g at 4°C for 1 h. The supernatant was 
collected and quantified with a 2-D Quant kit (GE Health- 
care). 

2D gel experiments and data analysis 

A series of 2D gels were produced essentially as reported 
previously (31). Isoelectric focusing was performed with the 
Ettan IPGphor 3 Isoelectric Focusing System. Immobiline 
nonlinear pH 3-10 DryStrips (GE Healthcare) were run us- 
ing rehydration buffer (8-M urea, 2% CHAPS and 20-mM 
Dithiothreitol, DTT) containing 0.5% v/v IPG Buffer (GE 
Healthcare). Sodium dodecyl sulphate-polyacrylamide gel 
electrophoresis (SDS-PAGE) was performed using 12.5% 
polyacrylamide gels without a stacking gel in the Ettan Dalt 
Six Elect Unit 230 (GE Healthcare). Gels were stained with 
0.04% w/v PhastGel Blue R (CBB R-350; GE Healthcare) 
in 10% acetic acid. In order to obtain reliable results from 
2D images, the experiments were performed independently 
four times. The protein spots were excised from 2D gels and 
subjected to tryptic digestion and liquid chromatography- 
tandem mass spectrometry (LC-MS/MS). 

Generation of antiserum against MET-8 

GST-MET-8 (containing MET-8 amino acids 575-689) fu- 
sion protein was expressed in BL21 cells, and the soluble 
recombinant protein was purified and used as the antigen 
to generate rabbit polyclonal antiserum, as described previ- 
ously (17,18). 

Protein analyses 

Protein extraction, quantification, western blot analysis and 
protein degradation assays were performed as described 
previously (32). Equal amounts of total protein (10 [xg) were 
loaded into each protein lane. After electrophoresis, pro- 
teins were transferred onto PVDF membrane. Western blot 
analyses were performed by using antibodies against the 
proteins of interest. 

Northern blot analysis 

Ribonucleic acid (RNA) was extracted as described previ- 
ously (33), and equal amounts of total RNA (20 |xg) were 
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loaded into each lane of 1 .3% agarose gels for electrophore- 
sis. The RNAs were transferred from gels to nylon mem- 
brane. The membranes were probed with RNA probes spe- 
cific to met-8, NCU06171, NCU0651 1 and mig-6 mRNA. 

Chromatin immunoprecipitation analysis 

Chromatin immunoprecipitation (ChIP) assays were per- 
formed as described previously (17). Briefly, crassa tis- 
sues harvested under the experimental conditions were fixed 
with 1% formaldehyde for 15 min at 25° C with shaking and 
then stopped with glycine at a final concentration of 125 
mM. Cross-linked tissues were ground and resuspended at 1 
g/ml in lysis buffer containing 1-mM PMSF, l-(jLg/ml pep- 
statin A and l-|xg/ml leupeptin. Chromatin was sheared 
by sonication to ~500-1000-bp fragments. One milliliter of 
protein (2-mg protein/ml) was used per immunoprecipita- 
tion, and 10 |xl was kept as the input DNA. ChIP was car- 
ried out with 10 |jl1 of antibody to trimethylated H3K9 (07- 
442; Millipore), 5 |jl1 of antibody to H3 (2650; CST), 2 |xl 
of antibody to trimethylated H3K4 (07^73; Millipore) or 
10 |jl1 of antibody 8WG16 to Pol II carboxyl-terminal do- 
main (CTD) repeats (ab817; Abeam). Immunoprecipitated 
DNA was quantified using real-time polymerase chain reac- 
tion (PCR) (CFX96; Bio-Rad) with primer pairs (see Sup- 
plementary Table SI). ChlP-quantitative PCR (qPCR) data 
were normalized by the input and presented as a percent- 
age of input DNA. Each experiment was independently per- 
formed at least three times. 

Western blot analyses of histone H3 and trimethylated H3K9 

Approximately 0.1 g of dried frozen tissue was resuspended 
in 1 ml of extraction buffer containing protease inhibitors. 
The mixture was vortexed and sonicated briefly. The ho- 
mogenate was then aliquoted (100 fxl) and centrifuged at 
10 000 g for 15 min at 4°C. The precipitate was resus- 
pended in 100 |jl1 of loading buffer (50-mM Tris pH 6.8, 
2% SDS, 1-mM ethylenediaminetetraacetic acid, 1% glyc- 
erol, 3-M urea, 100-mM DTT and 0.25% bromophenol 
blue) and boiled for 30 min. The crude extracts were cen- 
trifuged at 10 000 g for 10 min, and the supernatants were 
loaded onto 12% polyacrylamide gels for electrophoresis. 
After electrophoresis, proteins were transferred onto PVDF 
membranes and western blot analysis was performed using 
antibodies against trimethylated H3K9 (07^42; Millipore) 
and H3 (2650, CST), respectively. 

RESULTS 

MET-8 protein expression is decreased in strains with dele- 
tions in the Cul4-based histone methyltransferase complex 

Previous studies showed that the DCDC (the DIM-5/-7/- 
9, CUL4/DDB1 complex) (19) or Cul4-DDB1-DCAF26- 
DIM-7-DIM-5 complex is essential for histone H3K9me3 
and DNA methylation in N crassa (17,18). To investi- 
gate how the complex affects gene expression, we examined 
changes in the proteome of cul4^^ and ddbl^^ strains dur- 
ing vegetative growth in liquid medium under constant hght 
conditions. After 2D gel electrophoresis of soluble protein 
extracts, LC-MS/MS analysis of excised gel spots led to 



the identification of 172 spots representing 170 unique pro- 
teins with differences in intensity in the cul4^^ and ddbl^^ 
strains compared to the wild-type strain (Supplementary 
Figure SI). Among the spots that differed in intensity be- 
tween wild-type and mutant strains, MET-8 (NCU06512) 
was quite prominent, because its expression was signifi- 
cantly decreased in the mutant strains (Figure lA and B). 
MET-8 is responsible for transfer of the methyl group of 
5-methyltetrahydrofolate to homocysteine, resulting in me- 
thionine for the SAM cycle in N crassa (34). Three pro- 
tein spots, labeled in Figure lA and B, were identified by 
LC-MS/MS as the same MET-8 protein (Figure IC). To 
further confirm that these three protein spots were MET-8, 
we generated a MET-8 polyclonal antibody and performed 
western blotting to probe the PVDF membrane that held 
proteins transferred from the 2D gel of the wild-type sam- 
ple. As shown in Figure ID, the MET-8 antibody signal fits 
with the three protein spots identified as MET-8 by LC- 
MS/MS, which suggests that the protein spots may repre- 
sent post-translational modifications or other isoforms of 
MET-8. 

Since the Cul4-based histone methyltransferase complex 
is essential for H3K9me3 and heterochromatin formation 
in N crassa, we next tested whether the other members of 
this complex were involved in regulation of MET-8. To com- 
pare levels of MET-8 proteins among the wild-type, cul4^^ 
and ddbl^^ strains, we performed side-by-side comparisons 
of cultures harvested after growing under the same condi- 
tions. As shown in Figure 2 A, protein levels of MET-8 were 
lower in each of the mutants compared to the wild-type 
strain, indicating that the Cul4-based histone methyltrans- 
ferase complex was a positive regulator of MET-8 expres- 
sion. The low MET-8 levels in the mutants suggest that ei- 
ther MET-8 protein was less stable or transcription of met-8 
was decreased in the absence of each subunit. To test these 
possibilities, we compared the stability of MET-8 in wild- 
type, cul4^^ and ddbl^^ strains after addition of cyclohex- 
imide (CHX). MET-8 was very stable in the wild-type and 
mutant strains (Figure 2B and C), indicating that the low 
levels of MET-8 in mutants were not due to an increase in 
the degradation of MET-8 protein. We next measured met-8 
mRNA in the mutant and wild-type strains by northern blot 
analysis. As shown in Figure 2D and Supplementary Fig- 
ure S2, northern blot analyses revealed that levels of met-8 
mRNA were reduced in each of the deletion strains com- 
pared to the wild-type strain, indicating that the expression 
of met-8 was affected in these mutants. 

Since the Cul4-based histone methyltransferase complex 
generally plays an important role in regulating histone 
H3K9 trimethylation at AT-rich DNA regions that are po- 
tentially packaged into heterochromatin (5,6,17-19), these 
results suggest that defects in heterochromatin formation 
may affect the proper expression of met-8. 

Proper formation of heterochromatin is required for high ex- 
pression of met-8 

A crucial step in heterochromatin assembly is the binding of 
HPl to the trimethylated H3K9 produced by DIM-5. Sev- 
eral studies have shown that HPl proteins play important 
roles in heterochromatin formation and regulation of gene 
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Figure 1. Expression of methionine synthase (MET-8) is decreased in cul4^^ and ddbl^^ strains. (A) Total protein extracts from N. crassa wild-type (WT) 
cells were separated by 2D gel electrophoresis. The rectangle shows the location of MET-8 proteins in the 2D map. (B) Expression pattern of MET-8 proteins 
in WT, cul4^^ and ddbl^^ strains. The arrows indicate the spots of MET-8 proteins identified by LC-MS/MS. (C) Identification of MET-8 proteins in N. 
crassa using LS-MS/MS. E: experimental values obtained directly from gels; T: theoretical values obtained from ExPASy search for full-length precursor 
proteins; MW: molecular weight (kDa). (D) Protein spots of MET-8 in 2D gels were confirmed by western blot analysis using antiserum against MET-8. 
Five protein spots were found on a PVDF membrane (top panel) containing proteins transferred from a 2D gel with WT sample (bottom panel) stained 
by CBB R-350. 



expression, but they contribute little to H3K9me3 (14,35). 
We compared levels of MET-8 between the wild-type strain 
and the hpo deletion mutant, which is deficient for HPl. 
MET-8 levels were decreased in the hpo^^ strain compared 
to the wild-type and ku70^^^ strains (Figure 3 A). To further 
investigate whether proper formation of heterochromatin is 
required for high expression of met-8 in N crassa, we mea- 
sured the levels of met-8 mRNA in the hpo^^, wild-type, 
and ku70^^ stains. As in the cul4^^ strain, levels of met- 
8 mRNA in the hpo^^ strain were dramatically decreased 
compared to the wild-type and ku70^^ strains (Figure 3B 
and Supplementary Figure S2). In crassa, HPl is known 
to be involved in DNA methylation by recruiting the DNA 
methyltransferase DIM-2 to chromatin (14). To examine 
whether proper DNA methylation is required for high ex- 
pression of met-8 in N crassa, we measured the levels of 
MET-8 protein and met-8 mRNA in dim-2^^ and wild-type 
strains. Both MET-8 protein levels (Figure 3C) and met-8 
mRNA levels (Figure 3D) in the dim-2^^ strain were com- 
parable to those in the wild-type strain. Together, these re- 
sults indicate that proper formation of heterochromatin, not 



proper DNA methylation, is required for expression of met- 
8'mN crassa. 



Proper establishment of H3K9me3 is required for the highly 
efficient expression of met-8 

Previous studies showed that N. crassa DIM-5 is essen- 
tial for histone H3K9me3 and heterochromatin formation 
(15,19,35). To determine whether this modification was in- 
volved in maintaining high levels of met-8 expression, we 
generated a mutant strain with an amino acid substitution 
(lysine 9 to arginine; hereafter referred to as the H3K9R 
mutant). This mutation mimics loss of trimethylation on 
H3K9. The N. crassa genome has a single histone H3 gene 
(NCU01635) (36), which facilitates testing of mutations in 
this gene. To obtain an H3K9R mutant, we made a knock- 
in cassette containing the mutated histone H3 gene and a 
hygromycin resistance gene {hph) inserted downstream of 
the H3 gene J untranslated region (UTR) (Figure 4A). 
This cassette was then transformed into a ku70^^^ genetic 
background strain, and /z/?/z-resistant transformants were 
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Figure 3. Proper formation of heterochromatin is required for the highest 
expression of met-8. (A) Results of western blot analysis of MET-8 (top 
panel) and the densitometric analyses from four independent experiments 
(bottom panel) in WT, ku70^^^ and hpo^^ strains. (B) Northern blot anal- 
ysis showing met-8 mRNA in WT, ku70^^^ and hpo^^ strains. (C) Results 
of western blot analysis of MET-8 (top panel) and the densitometric analy- 
ses from four independent experiments (bottom panel). (D) Northern blot 
analysis showing met-8 mRNA in WT and dim-2^^ strains. 



selected on a plate with hygromycin B. PCR analysis con- 
firmed the integration of the knock-in cassette at the en- 
dogenous H3 gene locus, and the homokaryotic nature of 
the H3K9R strains was ensured by microconidia purifica- 
tion and germination of ascospores from crossing H3K9R 
mutant {bd, a) to 301-6 strain (bd, A) (Supplementary Fig- 
ure S3). DNA sequencing of the endogenous H3 gene fur- 
ther confirmed the mutated H3K9R and the presence of 
a nonsense mutation for the introduction of a Sa/I site in 
these strains (Figure 4B and Supplementary Figure S3A). 
We tested the levels of histone H3 or trimethylated H3K9 
using western blot analysis in wild-type, ku70^^, cul4^^, 
dim-5^^ and H3K9R strains. As shown in Figure 4C, all 
of the strains had similar levels of histone H3. In contrast 
to the robust H3K9 trimethylation in the wild-type and 
ku70^^^ strains, no trimethylated H3K9 was detected in the 
H3K9R mutant or in the cul4^^ and dim-5^^ strains (Fig- 
ure 4C). This confirmed that the H3K9R mutant strains had 
lost H3K9me3. Moreover, the H3K9R mutants exhibited 
dense, cauliflower-like growth patterns with abnormal hy- 
phae and asexual spores on plates, similar to the phenotypes 
of the cul4^^ and dim-5^^ strains (Figure 4D and Supple- 
mentary Figure S3B) (17,18). The H3K9R strain exhibited 
slow growth rates on race tubes compared to the wild-type 
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strain (Figure 4E), the growth rate being similar to that of 
the cul4^^ and dim-5^^ strains. These resuhs suggest that 
the lost H3K9 trimethylation contributed to most of the de- 
fective phenotypes we observed in strains with mutations 
in the Cul4-based histone methyltransferase complex in 
crassa. Next, we measured levels of met-8 mRNA in the 
H3K9R, ku70^^ and wild-type strains. Northern blot anal- 
yses showed that the level of met-8 mRNA in the H3K9R 
mutants was dramatically decreased compared to those in 
the ku70^^^ and wild-type strains (Figure 4F and Supple- 
mentary Figure S3C), demonstrating that H3K9 trimethy- 
lation is required for the high expression of met-8 in N 
crassa. 



Location and organization of the met-8 gene in linkage group 
III of A^. crassa 

To determine why the proper formation of heterochromatin 
is required for the high-efficiency expression of met-8 in N. 
crassa, we explored its chromosomal location and nearby 
DNA sequence features in the N. crassa database (http:// 
www.broadinstitute.org/annotation/genome/neurospora/). 
We found that the met-8 (NCU06512) gene was located 
in the left arm of linkage group III (located in Super- 
contig 3: 1698305-1701771) between the NCU06511 and 
NCU06171 genes. Unexpectedly, immediately upstream of 
the met-8 promoter, there was an ~50-kb AT-rich DNA 
region located between met-8 and NCU06171 (Figure 5 A 
and B). In crassa, histone H3K9 is always trimethylated 
at AT-rich DNA regions (35), thus histone H3K9 at this 
50-kb AT-rich region upstream of met-8 is likely to be 
trimethylated. We carried out a series of ChIP assays using 
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Figure 5. Localization and organization of the met-8 gene on linkage group III of TV; crassa. (A) and (B) Schematic depicting a 50-kb AT-rich DNA region 
located between met-8 and NCU06171 on linkage group III of crassa genome. Bars (primer pairs 1-5) above the schematic indicate the regions tested 
by qPCR in ChIP assays. (C)-(G) ChIP analyses showing enrichment of H3K9 trimethylation at different regions upstream oi met-8 in wild-type, cul4^^, 
ddbl^^, dcaf26^^, dim-7^^, dim-5^^ and hpo^^ strains. The loci detected by ChIP are indicated in (B). (H) and (I) Northern blot analysis showing levels 
of NCU06171 and NCU06511 mRNA in wild-type, kulO^^^ , cuM^^, ddbl^^, dcaf26^o, dim-?^^, dim-S^^ ^^^ko strains. 



antisera against trimethylated H3K9 and oligonucleotide 
primer pairs (1-5; see Figure 5B) designed to target the 
area from the met-8 promoter region to the middle of the 
proximal 50-kb AT-rich region in the wild-type strain and 
each mutant strain. qPCR with primer pairs 1 or 2 and 
immunoprecipitated DNA from the wild-type or any of the 
mutant strains did not generate product signals, indicating 
that there was no H3K9me3 at the met-8 promoter region 
(Figure 5C and D). However, as revealed by primer pairs 
3-5, H3K9me3 was dramatically enriched at the 50-kb 
AT-rich region in the wild-type and hpo^^ strains, but not 
in the other mutants (Figure 5E-G). This result indicates 



that histone H3K9 trimethylation of the 50-kb AT-rich 
region began at the region between primer pairs 2 and 3. 
Our Chip data demonstrate that the 50-kb AT-rich DNA 
region was modified by histone H3K9 trimethylation and 
that this region may form heterochromatin. 

To confirm whether proper heterochromatin formation in 
the 50-kb AT-rich region is required for high met-8 expres- 
sion specifically, we measured the expression of NCU0651 1 
and NCU06171, which are located in the vicinity of the 
50-kb AT-rich region. Northern blot analyses showed that 
the levels of NCU06511 and NCU06171 mRNA expressed 
in the mutants were comparable to those in the wild-type 
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Figure 6. High expression of the mig-6 gene depends on its proximal 4.2- 
kb AT-rich heterochromatin region. (A) Schematic depicting a 4.2-kb AT- 
rich DNA region (peak3) located between mig-6 and NCU 12057 on link- 
age group VII of TV crassa genome. (B) Northern blot analysis showing lev- 
els mig-6 mRNA in WT, kulO^^^ , cul4^^ , ddbl^^, dcaf26^^, dim-7^^ 
and dim-5^^ strains. (C) Northern blot analysis showing levels of mig-6 
mRNA in the WT and AT-rich region peak3 deletion strains. 



strain (Figure 5H and I), indicating that proper formation 
of the 50-kb AT-rich heterochromatin region only affects 
met-8 expression. Taken together, these results strongly sug- 
gest that the 50-kb AT-rich heterochromatin region may 
play an important role in regulating met-8 expression. 



Heterochromatin positively regulates mig-6 gene expression 

We tried to delete the 50-kb AT-rich DNA region at the 
endogenous location by replacement with the hph gene to 
see whether it was indeed responsible for activating met-8 
expression. However, we could only obtain heterokaryon 
strains, which showed no visible change in met-8 mRNA 
levels. Next, we wondered whether other heterochromatins 
can activate proximal gene expression other than the met- 
8 gene in the Neurospora genome. After screening a set 
of heterochromatin-nearby genes, we found the mig-6 gene 
(NCU09285) located nearby a 4.2-kb AT-rich region on 
linkage group VII (Figure 6A). A previous study showed 
that this 4.2-kb AT-rich region (referred to as peak3) en- 
riches peaks of H3K9me3 and DNA methylation (35). 
Northern blot analyses revealed that levels of mig-6 mRNA 
were reduced in each of the heterochromatin-defective 
strains compared to wild-type (Figure 6B). Moreover, dele- 
tion of its upstream heterochromatin region led to a de- 
crease of mig-6 mRNA levels compared to that in the wild- 
type strain (Figure 6C). Consistent with the regulation of 
met-8 expression, these results confirm that the 4.2-kb AT- 
rich heterochromatin region positively regulates mig-6 gene 
expression. 



Production of ectopically expressed met-8 is inefficient at the 
his-3 locus but efficient at downstream of heterochromatic 
regions 

To investigate a direct correlation between met-8 expression 
and its proximal 50-kb AT-rich heterochromatin region, we 
checked met-8 expression in an environment lacking this re- 
gion by cloning the entire met-8 gene, integrating it into 
the his-3 locus and detecting the output in different mu- 
tant strains. Bioinformatics analysis of the met-8 promoter 
by TESS (37) showed that its core region was about 1.0 
kb from the start codon to the boundary (between primer 
pairs 2 and 3) of the proximal 50-kb AT-rich region (Fig- 
ure 7A). To confirm this prediction, ChIP assays with an- 
tibody 8WG16 against the CTD repeats of the Pol II Rbpl 
subunit were performed to measure the enrichment of Pol 
II on the met-8 promoter. As shown in Figure 7B, enrich- 
ment of Pol II was dramatically increased at the primer pair 
1 region, but not at the other regions. Consequently, the 
1.0-kb core region and the 1.0-kb core region plus 500 bp 
of AT-rich sequences were cloned as promoters for ectopic 
met-8 expression. Full-length open reading frames and the 
y UTR of met-8 were amplified from genomic DNA and 
cloned together into both the pmet-8(1.0-kb)-5Myc-6His 
and pmet-8(1.5-kb)-5Myc-6His vectors (Figure 7C). We 
then introduced pmet-8(1.0-kb)- or pmet-8(1.5-kb)-driven 
Myc-tagged MET-8 constructs into the his-3 locus of 301- 
6 (bd, his-3, A) and dim-7^^ (bd, his-3) strains. As shown 
in Figure 7D, the expression levels of Myc-MET-8 at the 
his-3 locus of wild-type or dim-7^^ strains driven by met-8 
promoters of different lengths were comparable, indicating 
that the defect in heterochromatin formation had no effect 
on ectopic expression of met-8. Next, we performed west- 
ern blot analyses using antiserum against MET-8 proteins 
in these transformants to compare expression between the 
ectopic met-8 gene and endogenous met-8 gene in wild-type 
and dim-7^^ strains. As shown in Figure 7E, the outputs of 
ectopic MET-8 were much lower than endogenous MET- 
8. Densitometry analyses revealed that the difference in 
protein expression between ectopic and endogenous MET- 
8 was at least 3 -fold in the wild-type background (Figure 
7F). Interestingly, when we introduced the pmet-8(1.0-kb)- 
driven Myc-tagged MET-8 construct with a 4-kb AT-rich 
sequence from its upstream 50-kb AT-rich region into the 
his-3 locus of 301-6 (bd, his-3. A) strain, the outputs of 
ectopic myc-met-8 gene in transformants were higher than 
those of the pmet-8(1.0-kb)-driven Myc-MET-8 construct 
(Supplementary Figure S4), suggesting that the ectopic ex- 
pression of met-8 genes was positively regulated by this 4-kb 
AT-rich region. 

We also wondered whether the expression of the met-8 
gene can be stimulated by other heterochromatin regions. 
To test this possibility, we examined the Neurospora genome 
and focused on two AT-rich DNA regions located in linkage 
group III (Figure 8A, top panel) and linkage group I (Fig- 
ure 8A, bottom panel), respectively. ChIP assays confirmed 
the enrichment of trimethylated H3K9 in these regions (Fig- 
ure 8B), indicating that they are potential heterochromatin 
domains. 

We next generated two knock-in cassettes contain- 
ing pmet-8(1.0-kb)-5Myc-6His-MET-8 and the selectable 
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marker hph targeting downstream of each previously men- 
tioned region. These cassettes were then transformed 
into a ku70^^^ genetic background strain, respectively. 
The homokaryon pmet-8(1.0-kb)-5Myc-6His-MET-8::/z/?/z 
strains integrated downstream of AT-rich region 1 or re- 
gion 2 were obtained by microconidia purification (here- 
after referred to as ATl-MET-8 and AT2-MET-8, respec- 
tively) (Figure 8C). As shown in Figure 8D and E, the 
outputs of ectopic met-8 gene downstream of AT-rich re- 
gion 1 or region 2 were much higher than those at the 
his-3 locus, indicating that the ectopic expression of met- 
8 genes was positively regulated by these two proximal 



heterochromatin domains. Interestingly, the outputs of ec- 
topic myc-met-8 gene downstream of AT-rich region 1 or 
region 2 in ddbl^^ strains were lower than those in ku70^^^ 
genetic background strains (Supplementary Figure S5), fur- 
ther confirming that the ectopic expression of met-8 genes 
was positively regulated by the proper formation of its prox- 
imal heterochromatin domains. 

Taken together, our experiments indicate that peak pro- 
duction of protein encoded by the met-8 gene is dependent 
on its proximal AT-rich heterochromatin regions. 
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Defects in heterochromatin formation affect H3K4me3 and 
recruitment of RNA Pol II at the met-8 5' UTR 

To provide a possible mechanistic explanation for how met- 
8 expression is regulated by its proximal 50-kb AT-rich het- 
erochromatin, we first checked the density of histone H3 at 
the met-8 promoter region in wild-type and mutant strains. 
Chip assays showed that the density of H3 at the met-8 
promoter region was comparable in wild-type and mutant 
strains (Figure 9A), indicating that the defect in heterochro- 
matin formation had no effect on nucleosome assembly at 
the met-8 5' UTR. 

Histone H3 lysine 4 trimethylation (H3K4me3) occurs 
in the promoter regions and 5^ UTR region of actively 
transcribed genes, and it is associated with transcription 
initiation. We next carried out ChIP assays to compare 
H3K4me3 at the met-8 5^ UTR region in wild-type and mu- 
tant strains. As shown in Figure 9B, levels of H3K4me3 



were enriched at the 5^ UTR region of met-8 in the wild-type 
strain, whereas levels of this modification were decreased at 
this region in mutant strains, indicating that histone modifi- 
cations at met-8 were altered in mutant strains. These results 
suggest that the defect in met-8 proximal heterochromatin 
affected histone modifications of the met-8 gene. Consider- 
ing the distribution of H3K4me3, which serves as a marker 
of active genes, these results may explain the downregula- 
tion of met-8 in mutants. ChIP assays with antibody 8WG16 
were carried out in wild-type and mutant strains to examine 
the enrichment of RNA Pol II at the 5' UTR of met-8. In the 
absence of Cul4 or each of the other subunits, enrichment of 
RNA Pol II at met-8 promoter regions was lower compared 
to wild-type strains (Figure 9C), suggesting that the defect 
in heterochromatin formation affected the recruitment of 
RNA Pol II to the met-8 promoter in mutants. Consistent 
with the decreased expression of met-8 in the mutant strains, 
these results indicate that proper formation of heterochro- 
matin is required for proper histone modifications of met-8 
and recruitment of RNA Pol II to the met-8 promoter. 

DISCUSSION 

N. crassa, a filamentous fungus, is a powerful, genetically 
tractable system for studies of heterochromatin formation 
and gene expression. In this study, we found that a het- 
erochromatic region was necessary for the high-efficiency 
expression of met-8, contradicting the well-accepted no- 
tion that heterochromatic regions create repressive environ- 
ments for gene silencing by recruiting certain components 
(38). The met-8 gene has an ~ 50-kb AT-rich DNA region 
contiguous with its promoter, which forms heterochromatin 
and positively regulates met-8 expression. In the absence of 
the proximal 50-kb AT-rich DNA region, expression of met- 
8 was dramatically decreased. Furthermore, we showed that 
mig-6 expression was directly regulated by its proximal 4.2- 
kb heterochromatin region. To the best of our knowledge, 
this is the first systematic functional analysis of heterochro- 
matin positively regulating proximal gene expression in fil- 
amentous fungi. 

In D. melanogaster, Howe et al. demonstrated that the 
level of light gene expression was correlated with the size 
of the heterochromatin block adjacent to that gene (23). 
In fact, it is well known that heterochromatic genes in D. 
melanogaster are repressed when moved to distal euchro- 
matin by chromosomal rearrangements (39). This indicates 
that proximity to heterochromatin is an important regula- 
tory requirement for the normal functioning of heterochro- 
matic genes (22,25). Consistent with the observations in D. 
melanogaster, our results demonstrate that expression of the 
met-8 or mig-6 gene relies on their proximal heterochro- 
matic region in crassa. 

In D. melanogaster, the expression of several well-studied 
heterochromatic genes, including light and rolled, is com- 
promised in a genetic background in which the HPl (named 
SU(VAR)205 in D. melanogaster) has been reduced (26,40). 
Here, we show that peak expression of met-8 was compro- 
mised in a genetic background deficient for HPl (hpo^^ 
strain) in crassa. A major challenge is to understand 
how the molecular functions of SU(VAR) proteins that 
affect heterochromatic gene expression can be reconciled 



6194 Nucleic Acids Research, 2014, Vol 42, No. 10 



with the current models of SU(VAR)- and heterochromatin- 
mediated euchromatic gene silencing. A large-scale map- 
ping analysis in Drosophila embryonic cells showed that 
HPl protein is distributed throughout the concertina, light 
and rolled heterochromatic gene regions and that it binds 
both unique and repetitive sequences in intronic and exonic 
portions of the genes, respectively (41,42). These data sug- 
gest that genes around heterochromatin may have evolved 
a transcriptional dependence on factors that are known to 
silence gene expression. However, most Su(var) mutations 
are pleiotropic and most SU(VAR) proteins are multifunc- 
tional, making it more difficult to determine a clear-cut ex- 
planation for these phenomena. 

To determine if the trimethylation at lysine 9 of H3 
plays an important role in met-8 transcription, we created a 
H3K9R substitution strain. As in the cul4^^ strain, expres- 
sion of met-8 was dramatically decreased in the H3K9R mu- 
tants, indicating that the proper formation of heterochro- 
matin is required for the highest-efficiency expression of 
met-8 in N crassa. Consistent with the idea that SU(VAR) 
proteins affect expression of these genes by direct binding, 
our results clearly show the effects of these modifiers and 
that H3K9me3 directly targets heterochromatic regions, 
which in turn influences met-8 expression. 

Another interesting question is how might the peak ex- 
pression of met-8 be compatible with the known silencing 
properties of heterochromatin, and what factors account 
for the difference between functional and silent heterochro- 
matin? 

As we have mentioned, methionine synthase is a key en- 
zyme for methionine synthesis in crassa, and highly ex- 
pressed met-8 is required throughout all phases of the cell 
cycle. One explanation, therefore, is that the requirement for 
the continuous activity of MET-8 and similar gene prod- 
ucts resulted in a robust adaptation strategy to promote 
continued expression in spite of repressive heterochromatic 
microenvironments. Our ChIP data revealed a sharp en- 
richment of H3K9me3 just beyond the met-8 promoter at 
the boundary of its proximal 50-kb AT-rich DNA region 
(between primers 2 and 3) in the wild-type strain, suggesting 
that there is a clear boundary between the met-8 promoter 
and its proximal heterochromatin. Therefore, we would ex- 
pect that strong boundaries are required both to maintain 
the integrity of the heterochromatin domains that facilitate 
heritable gene silencing and to protect neighboring gene loci 
from pervasive spreading of heterochromatin. 

In addition, in our study, although levels of ectopic ex- 
pression of met-8 were dramatically lower than those of 
the original met-8 gene, western blot analyses with c-Myc 
antibody revealed high levels of Myc-MET-8 under nor- 
mal growth conditions, indicating that the met-8 gene has a 
strong promoter, which is necessary to ensure that enough 
methionine synthase will be produced to generate suffi- 
cient methionine and SAM for cell viability. Several stud- 
ies in yeast, fruit fly and human cells have shown that ar- 
tificial recruitment of specific transcription activators, such 
as Rebl and Gal4, to DNA elements is sufficient to block 
the spread of heterochromatin (43). Moreover, strong RNA 
polymerase II promoters can also insulate genes in S. cere- 
visiae and S. pombe, again without requiring transcription 
(43,44). Therefore, another plausible explanation is that 



met-8 gene is insulated from repressive heterochromatin by 
its strong promoter, essentially maintaining a microenviron- 
ment that is compatible with transcription. 

Remarkably, our ChIP data revealed a sharp enrichment 
of RNA polymerase II at the met-8 promoter side of the 
met-8 /50-kb heterochromatin boundary, coincident with 
the increased enrichment of H3K4me3 at this region. Con- 
sistent with our observations, a truncated human CD2 gene 
lacking part of its locus control region (LCR) showed varie- 
gated expression in the thymus when the transgene was in- 
tegrated into or close to pericentromeric heterochromatin 
in mice (45). Surprisingly, a CD2 gene that contained the 
complete CD2 LCR yielded full nonvariegated expression 
when integrated into heterochromatin. Although this is a 
particularly striking result, since it contradicts the long-held 
view that heterochromatin has a dominant silencing effect 
on euchromatic genes that become associated with it, these 
results strongly indicate that insulation at these sites is a 
consequence of direct competition between the spreading 
heterochromatin and transcription factor binding coupled 
with chromatin remodeling at the barrier (46). 

CONCLUSION 

The findings of heterochromatin-mediated stimulation of 
gene expression expand our understanding of the function 
of heterochromatin and demonstrate that heterochromatic 
genes are evolutionarily conserved in eukaryotes. They also 
emphasize the importance of the chromosomal context in 
complex eukaryotic genomes. We predict that the mech- 
anisms of heterochromatin-mediated stimulation of gene 
expression in eukaryotes will soon be revealed in future 
studies. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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